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APPLICATION OF HUMAN TRANSFER FUNCTIONS 


TO SYSTEM ANALYSIS 

By James J. Adams and Maxwell W. Goode 
Langley Research Center 

SUMMARY 

An analytical study was made of a full-scale, manually controlled, lunar-landing 
simulator by using analytical transfer functions for the pilot control response along with 
the analytical representation for the mechanisms. The simulator reproduced lunar grav- 
ity by supporting five-sixths of the weight of the test vehicle with an overhead cable. The 
cable was kept directly over the test vehicle by the automatic control of the longitudinal 
drive mechanism of the simulator. The results showed that some of the dynamic charac- 
teristics of the simulator were in a range to influence the response of the manually con- 
trolled systems which were to be tested. 

The results of the analytical study were checked by actual operation of the simula- 
tor. The simulator was operated with the gain of the longitudinal drive set as high as was 
possible with the actual mechanism and with a lower gain to determine whether this change 
would affect the pilot's response. The pilots reported that the slower responding system 
(low gain) was more difficult to control. The records clearly showed a decrease in 
system damping with the degraded system. 

INTRODUCTION 

One of the reasons for determining human transfer functions is to allow evaluation 
and prediction of the performance of manually controlled systems during design studies. 
This report explains the use of such transfer functions in an analysis of the drive system 
of a lunar-landing simulator to determine the effect of the drive system characteristics 
on the performance of the pilot -vehicle combination. 

The lunar -landing simulator is a manually controlled six-degree -of -freedom 
moving-base, lunar-gravity simulation device designed to define desirable handling quali- 
ties for a lunar-landing vehicle and piloting problems associated with lunar-landing 
maneuvers. Since no direct comparison could be made between the simulator and the 
actual cases, it was desirable to define the influence of the simulator on the piloting tasks. 
In this effort analytical transfer functions for pilot control responses and an analytical 
representation of the drive mechanisms were used. 
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The transfer functions used to describe the pilot's control action were derived in 
reference 1. These transfer functions describe the control used by the pilot when con- 
trolling a multi-loop system representative of the lunar -landing horizontal-translation 
maneuver. The use of these pilot transfer functions in determining the most suitable 
simulator drive characteristics is presented. 
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SYMBOLS 


equivalent leakage orifice area, ff2 (m2) 

relief valve area, ft^ (m2) 

speed of propagation, ft/sec (m/sec) 

in-lb ( m-N 


load viscous friction, 


rad/sec \rad/sec 


motor damping, — - n ^ 

rad/sec \rad/sec 

motor displacement, in3 (cm3) 
pump displacement, in^ (cm3) 


distance from whiffletree to vehicle, ft (m) 
gravity, 32.2 ft/sec2 (9.81 m/sec2) 


maximum deflection of first cable vibration mode, ft (m) 
load inertia, in-lb-sec2 (m-N-sec2) 
motor inertia, in-lb-sec2 (m-N-sec2) 
leakage coefficient, 

psi \ N/ m2 j 

torque coefficient, - - ( ■ m ~N \ 

rad/ sec \rad/ sec/ 

orifice flow coefficient, in 2 -lb 1 / 2 / m? . -N 1 / 2 ] 

’ sec \ sec ] 

tire torsional spring constant, in-lb/rad (m-N/rad) 
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V K <? 
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Pm 

Qc 

Ql 

Qm 

Qp 

Qrl 

R 
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T 

T L 

T m 
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V 


gains in analytical transfer function of pilot 
simulator drive system control gains 
constant 

pendulum length, ft (m) 
whiffle tree coordinate, ft (m) 
mass, slugs (kg) 
gear ratio 

motor pressure, lb/in2 (N/m2) 
compressible flow, in^/sec (cm^/sec) 
leakage flow, in3/sec (cm^/sec) 
motor flow, in^/sec (cm^/sec) 
pump flow, in*Vsec (cm^/sec) 
relief valve flow, in^/sec (cm^/sec) 
tire radius, in. (cm) 

Laplace operator, per second 
tension, lb (N) 
tire torque, in -lb (N-m) 
induction motor torque, in-lb (m-N) 
time, sec 

oil volume under compression, in^ (cm3) 
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cable density, lb/ft (N/m) 
translation, ft (m) 
cable coordinates, ft (m) 
bulk modulus of oil, lb/in2 (N/m^) 

control moment, rad/sec^ 

pitch-attitude angle, radians 

filter breakpoint frequency, rad/sec 

pendulum angle, radians 

motor angle, rad 

load shaft angle, rad 

motor speed, rad/sec 


0 ) n undamped natural frequency, rad/sec 

o)p pump speed, rad/sec 

k’sync synchronous speed, rad/sec 

Subscripts: 


B 

V 

c 

e 

f 

w 

m 


bridge 

vehicle 

command 

error 

filter 

whiffletree 

motor 
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Dots over symbols denote derivatives with respect to time. 


DESCRIPTION OF PROBLEM 


Simulator 

The description and operation of the lunar -landing simulator are given in detail in 
reference 2. The simulator itself is shown in figure 1. Basically, the simulator is a 
six-degree -of -freedom, moving base, lunar -gravity simulator with a flight envelope 
360 feet (110 m) by 42 feet (13 m) by 180 feet (55 m) high. Lunar gravity is simulated 
by supporting five-sixths of the weight of the test vehicle with a cable suspension system. 
The tension in the cables is regulated by measuring the load with a strain gage and oper- 
ating an overhead winch in response to the error in the measured tension. The suspen- 
sion cables are kept directly over the vehicle by measuring the angle of the cables with 
respect to the vertical and driving a traveling bridge in response to this measured angle. 
The longitudinal drive system of the traveling bridge is the subject of study in this 
report. See sketch (a). 
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Sketch (a).- Bridge and vehicle. 

Simplified equations defining the longitudinal motion of the overhead bridge and the 
cable -supported vehicle are: 

( m B + m VpB + m V l $ = 0 (!) 

myZXg + my l“(p + myg lep = 0 (2) 

The characteristic equation for the system in Laplace notation is 


m-o + mu 


rm 


ri 2 ) s 2 + ( m B + m v) ( m V gl ) 


m y 2 f 2 s 2 


= 0 


( 3 ) 


which is the familiar equation for a pendulum in which the frequency is primarily deter- 
mined by the length l but with an additional term representing the influence of the 
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pendulum pivot point being free to move and increasing the frequency somewhat. If bridge 
acceleration is a function of cp and cp, the equations of motion become 


(m B + niyjxg + my lip - K^cp + K^cp 

(4) 

my&B + m V^ 2< P + IHygZip = 0 

(5) 


the characteristic equation, in Laplace notation, being: 


(m B + my^myZ 2 js 2 + (fnyZK^Js + (m B + my^mygZj - my 2 Z 2 s 5 


m V lK <p 


= 0 ( 6 ) 


It can be seen from this equation that the gain, in the coefficient of s, will supply 

damping to the system, and the K ^ gain will increase the pendulum frequency and keep 
the bridge above the suspended vehicle. 


The bridge drive unit is an electrohydraulic unit consisting of a synchronous elec- 
tric motor which drives a variable -displacement hydraulic pump. Control of the bridge 
is exercised by the operation of a pump stroker, which controls the displacement of this 
pump. The pump drives fixed -displacement hydraulic motors attached to the wheels of 
the bridge. A fixed displacement of the stroker produces a steady-state constant veloc- 
ity of the bridge . 

Since the stroker controls bridge velocity instead of bridge acceleration, the control 
function of accelerating the bridge as a function of pendulum angle is achieved by dis- 
placing the stroker as a function of the integral of the pendulum angle, and the function of 
accelerating the bridge as a function of rate of change of pendulum angle is achieved by 
displacing the stroker as a function of pendulum angle. 


The predominant dynamic characteristic of the drive unit is the oscillatory response 
of bridge velocity to stroker displacement that results from the compressibility of the 
hydraulic fluid as it reacts against the mass of the bridge. This dynamic characteristic 
is expressed by the equation: 


X B 

Stroker displacement 


C l/ D m 


V J L ^2 


V 

13 D 


m 


+ ^s 

Dni 


+ 1 


(V 


where 

Cl 


0 


oil volume delivered by 
(cm 2 /unit) 

motor displacement, 28 
bulk modulus of oil, 1 x 


pump per unit displacement of stroker, 

in 2 /rad (460 cm 2 /rad) 

10 5 psi (6.9 x 10 8 N/ m 2 ) 


in 2 /unit 
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V 

Jl 


k l 


oil volume under compression, 600 in^ (9850 cm^) 

total inertia of load reflected at motor output shaft, 3050 in-lb-see^ 
(345 m-N-sec^) 

leakage coefficient in^/sec | cm3/sec\ 

P S1 \N/m 2 / 


For the system under study, the natural frequency of the drive unit as determined by 
this equation is 6.5 rad/sec. The detailed analysis of the system naturally included this 
characteristic. However, instead of using the transfer function presented in equation (7), 
a nonlinear representation for the hydraulics was used in the analysis. This representa- 
tion is presented in appendix A. 


Another important dynamic characteristic of the system is the oscillatory charac- 
teristic of the cable. Cable vibrations add to the measured cable angle and therefore add 
a spurious signal to the control signal. Equations for the first two modes of vibration for 
different fixed cable lengths were determined and used in the analysis. These equations 
included the effect of a lumped mass located near the vehicle representing the whiffletree 
which is a part of the vehicle suspension system. The derivation of these equations is 
given in appendix B. The frequency of the first mode for a 200-foot (61 m) cable length 
was 8.72 rad/sec which is very close to the natural frequency of the drive unit, and there- 
fore, it put a limit on the precision with which the bridge could track the vehicle. This 
vibration frequency is a function of cable length and vehicle weight. The analysis was 
made for fixed cable lengths and vehicle weights. 


Pilot Transfer Function 

The pilot transfer functions, which were used in conjunction with the lunar-landing 
simulator equations, are derived in reference 1 and are repeated here. These transfer 
functions were obtained from tests made with simplified fixed -base simulators. As is 
pointed out in reference 1, a combination of tests were used to obtain these transfer 
functions. Multi-axes compensatory tracking tests using automatic parameter tracking 
methods to determine the response of the pilot were used to obtain the inner-loop transfer 
functions. Multi-loop simulations and trial-and-error methods were used to obtain the 
outer-loop transfer functions. It should be noted that slightly different transfer functions 
have been measured with different pilots and that the particular transfer function selected 
for this study represents a nominal pilot response. 

The pilot-vehicle system involved in the landing maneuver is a multi-loop system 
described by the block diagram presented in figure 2. The inner loop deals with the atti- 
tude control of the vehicle. The vehicle response to attitude control was assumed to con- 
tain a proportional rate feedback and is given by the equation: 
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which defines a rate command system with a time constant of 2.0 seconds. Reference 1 
demonstrates that a pilot's response in such an inner loop is given by: 


KT fi + V\ 

6 _ 1 \ t / _ 96(1 + 0.4s) 

9e (s + t ) 2 (s + 6) 2 


(9) 


The combination of the pilot and vehicle gives this inner loop a closed-loop characteristic 
frequency of 1.2 rad/sec and a damping ratio of 0.26. 


The outer loop of the system deals with the longitudinal translation, and the vehicle 
response to attitude angle is given by a pure inertial response 


^V_ 5,36 

e s 2 


( 10 ) 


The relation is derived from the linearized equation of motion for the horizontal 
component of acceleration due to the l/6g thrust that would be in effect in the lunar 
environment, 

Xy = igsin0 (11) 

When small-angle linearization is used, 


Xy = ~ g8 = 5.360 ft/sec2 = 1.630 m/sec2 (12) 

_ k ft / sec2 _ i 63 m/sec2 

0 " ,3 rad ~ 1>63 rad (13) 

Reference 1 demonstrates that the pilot response in such an outer loop is 


0 C _ 0.9(1 + 9,2s) 
Xe (s + 10)2 


(14) 


When all these transfer functions are combined, block diagram algebra being used, a 
characteristic response of the complete system is obtained which has two small real 
roots, s = -0.167 and s = -0.336. In terms of an oscillatory response, these two roots 
define an overdamped response with a natural frequency given by 

u) n = V(0.167)(0.336) = 0.236 rad/sec (15) 

This system frequency characterizes the translation response of the system. Since the 
response characteristic of the longitudinal drive system of the simulator must have a 
response frequency higher than that of the system which is to be tested, this calculated 
response characteristic of the pilot-controlled translation response provides a first 
rough criterion for the required characteristic of the longitudinal drive system. 
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ANALYSIS AND TEST RESULTS 


Initial Study 

A detailed analytical study was made to determine the drive system characteristics 
and to evaluate the suitability of these characteristics for the purposes of the simulation. 
The complete computer program used in this study including the drive system, the cable 
dynamics, and the pilot -controlled vehicle is presented in figure 3. The system param- 
eters based on the physical constraints of the simulator were estimated. 

The first phase of the study was conducted to determine the response characteris- 
tics of the bridge. Open-loop step thrust inputs to the vehicle were used as the forcing 
function in these studies. The results showed that the presence of the cable vibration 
modes of motion in the system placed an upper limit on the pendulum damping gain K^,. 

If this gain was adjusted too high, the first vibration mode would become unstable, as is 
illustrated in figure 4. 

The limit on the pendulum damping gain placed further restriction on the pendulum 
frequency gain K<p. Representative system characteristics of vehicle velocity response 
to a 2 -second thrust impulse are shown in figure 5. The oscillatory nature of these 
responses is the result of the bridge drive system characteristics. It can be seen that a 
well -damped response with a frequency of 1.57 rad/sec (a period of 4 seconds) or a poorly 
damped response with a frequency of 2.5 rad/sec (a period of 2.5 seconds) could be 
achieved. 

Both of these frequencies are above the 0.236 rad/sec frequency for the pilot- 
controlled translation response of the lunar -landing system derived in the previous sec- 
tion. However, it cannot be confidently concluded that these frequencies are high enough 
to have no effect on the simulation. To determine the effect that the bridge response 
might have on the pilot-controlled maneuver, the analytical representation of the pilot and 
vehicle were included in a closed-loop representation of the complete system, and the 
response to a commanded 225-foot (69 m) displacement was determined. The results are 
presented in figure 6 which shows (1) the response of the pilot -vehicle combination alone, 
and (2) the response of the pilot-vehicle -simulator combination with the two different 
simulator characteristics presented before. By using the pilot -vehicle combination 
response as the standard for comparison, it can be seen that including the simulator 
bridge dynamics in the loop does indeed influence the response. With the lower gain 
(lower frequency) bridge control, the system is degraded to the point of instability. It 
was, therefore, concluded that the simulator should be operated with a frequency response 
as high as possible and that even then some influence of the simulator dynamics may be 
reflected in pilot control responses. 
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Post -Analysis Tests 


When the lunar -landing simulator was put into operation, the system characteris- 
tics which could be achieved with the actual mechanism were determined. It was found 
that the highest stable pendulum frequency that could be obtained was 1.4 rad/sec (a 
period of 4.5 seconds). By lowering the gain, a well-damped response with a fre- 

quency of 0.8 rad/sec (a period of 8 seconds) could be obtained. Time histories of cable 
angle response to open-loop thrust impulses with these two systems are shown in figure 7. 
The 4.5-second and 8-second period pendulum responses are clearly seen in this figure as 
well as the higher frequency cable vibration modes of motion. 

Runs in which a pilot closed the pitch and displacement loops were then made with 
both the high-gain and low -gain systems to check the results of the analytical study. 

These tests were started with the vehicle hovering at an altitude of approximately 30 feet 
(9 m). The pilot attempted to translate a prescribed distance and establish a hover over 
a mark located on the ground. Runs with the two systems were made one right after the 
other so that, insofar as possible, all environmental conditions other than the change in 
bridge response would be the same. 

It should be noted that the vehicle attitude -control system used in the simulator was 
different from the representation used in the analytical study. A rate feedback was incor- 
porated in the simulator, but the jet thrusters were operated in an on-off mode rather than 
in a proportional or modulated mode. The angular acceleration provided by the on-off 
system was 8°/sec 2 and the maximum angular velocity was 30°/sec. It is felt that the 
effects of the changes in the outer-loop bridge drive characteristics would be sufficiently 
well defined with each of the two attitude control systems so that the effects as obtained 
in the analytical study and in the simulator could be compared. 

Figure 8 is the recorded time history obtained in one of the lunar -landing simulator 
tests. In this test the first maneuver, from the 50-foot (15 m) mark to the 275-foot 
(84 m) mark, was done with the high-gain system. At the 50-second point, while the 
pilot was turning 180°, the drive system gain was adjusted to the lower setting, and the 
pilot started to go back toward the 50-foot (15 m) mark. The general nature of these 
maneuvers is very similar to that of the computed responses. With the high-gain system 
the pitch angle was well controlled, and with the low-gain system a decrease in the 
damping of the pitch-angle response can be seen at about 75 seconds. The pilot did not 
complete the translation in this maneuver, but landed when he stopped the pitch oscilla- 
tion because of a low fuel state. 

A series of tests was made at a later date than the tests presented in figure 8. 

Some additional refinements had been included in the bridge drive control system. Time 
histories of the cable -angle response to step commands in bridge acceleration for the 
refined bridge drive system are shown in figure 9. The periods of the bridge motion were 
5 and 8 seconds with the high and low gains, respectively. The damping characteristics 
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of the bridge were improved. Figure 10 is a reproduction of the time histories obtained 
in these tests in which a 100-foot (30 m) translation was performed, and additional time 
was spent in holding the hover condition. In this particular test the low -gain system was 
tested first, and then the high -gain system. Once again the decrease in the damping of 
the pitch-angle response can be clearly seen. 

All the pilots commented that they experienced an increase in difficulty in per- 
forming the maneuvering with the low -gain system. One pilot reported that he felt he 
was in ’’pilot -induced oscillation" condition with the low -gain system. 

Iterated Analytical Study 

Since the lunar-landing simulator characteristics that were achieved with the actual 
mechanism were not the same as those determined in the analytical study, and since the 
pilot's response showed a lower pitch -attitude angle limit in the flight tests than was 
assumed in the analytical study, the analytical study was repeated in an attempt to repro- 
duce more closely the flight time histories. In these repeated calculations the drive sys- 
tem gains were adjusted to give a poorly damped, 4.5-second-period pendulum mode of 
motion in one case, and a well-damped, 8-second period mode in the second case to corre- 
spond to the two conditions that were tested in flight. The same linear transfer functions 
were used for the representation of the pilot. The pitch -attitude limit, which is the limit 
placed at the output of the outer -loop analog pilot, was adjusted to a 10° value. This 10° 
limit corresponds more closely with the limit used by the pilots in the flight tests than 
did the 40° limit used in the initial analytical study. With this 10° pitch-angle limit, a 
225-foot (69 m) translation with a well-controlled attitude time history was obtained with 
the high-gain system, and this response is a better reproduction of the flight time history 
than was the initially calculated response. Once again a deterioration in damping, com- 
pared with that obtained with the high-gain system, was obtained with the low-gain sys- 
tem. These responses are presented in figure 11. 

There are some limitations to the pilot transfer functions as they are known at the 
present time. For example, in the particular case studied in this paper, a mistake was 
made in the assumed pitch-angle limit in the initial analytical study. Also, it is known 
that different pilots will control a given system with different transfer functions. The 
effects of these differences between subjects has not been studied in this paper. There 
may be other factors which are incorrect in the pilot transfer functions as they are known 
at the present time. However, useful information obtained in the example presented in 
this paper does add confidence to these functions. 

The general agreement between the analytical cases and the simulator cases as 
shown in figures 6, 8, 10, and 11 indicates that the human transfer functions are a rea- 
sonable representation of human response within limitations. It is felt that these transfer 
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functions are well enough defined to be employed advantageously in the design phases 
of similar manually controlled systems provided the limitations are recognized and 
are not exceeded. 


Single -Loop Task 

The analytical method, which was used to predict the effect of inserting a servo or 
servomechanism in the outer loop of a multi-loop problem in the previous sections of this 
paper, can also be used in single-loop or inner-loop control tasks. The primary differ- 
ences in the two problems are that single -loop control tasks are operated by human con- 
trollers at much higher frequencies, and that the response of the inserted servo must, 
correspondingly, be much higher in order to have no effect on the system. 

The analytically predicted effect of inserting first-order and second-order mecha- 
nisms in a single -axis control task are presented in figure 12(a). The vehicle dynamics 

n p 

considered here is a pure acceleration §- = The transfer function used to represent 

o s 2 

the pilot was taken from reference 1 and is 

5 _ 56(1 + 1.35s) 

0e (s + 7) 2 


The closed-loop response of this system contains an oscillatory characteristic with an 
undamped natural frequency of 3.3 rad/sec and a damping ratio of 0.36. 


The effect on the closed-loop characteristics of this system that results from adding 

T f 


a first-order lag to the dynamics of the form 


is also shown in figure 12(a). The 


S + Tf 

value of the breakpoint frequency tj was varied from 100 down to zero, and the change 
in system closed-loop damping ratio and natural frequency is plotted. Very drastic 
changes occur in the system response if the breakpoint frequency is reduced below 
20 rad/sec. 


The effect of adding a second-order servo to the dynamics is also shown. The form 

w n f 

of this inserted mechanism is s — . The assumed natural frequency 


s 2 + 2(0.7)co n f s + w n? f 


o> n was varied from 50 to 0, and the damping ratio was held constant at 0.7. The effect 
of inserting such a mechanism in the control loop is similar to that of the first-order lag, 
but more severe. Again, very drastic changes occur if the natural frequency is dropped 
below 20 rad/sec. 


Reference 3 presents results of an experimental study that is very similar to 
the single -loop analysis. In this experimental study the vehicle dynamics were repre- 
sentative of an aircraft longitudinal control system. A fixed-base simulation of a 
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landing -approach task was used for the control problem. A second-order representation 
of the elevator actuator was inserted in the control loop. The natural frequency of the 
actuator was varied from 62 rad/sec (10 cycles per second) down to 6.2 rad/sec (1 cycle 
per second). Pilot ratings were obtained for a number of natural frequencies within the 
range and are shown in figure 12(b). The damping of the second-order response was 
held constant at a damping ratio of 0.7. Very drastic changes in pilot ratings were 
obtained as the natural frequency was reduced below the 18 to 30 rad/sec range. Typical 
ratings were 2.5 (good, acceptable) for natural frequencies of 18 rad/sec and above, and 
6 (unsatisfactory) at a natural frequency of 6 rad/sec. From these results it can be seen 
that the experimental study confirms the analytical study of the single -loop, or inner - 
loop, problem. 


CONCLUDING REMARKS 

The results of this study have indicated that the human transfer functions, subject 
to the conditions of this experiment, are a reasonable representation of human response. 
These transfer functions appear to be defined well enough to be used to establish system 
specifications in the design phases of systems similar to those in the two examples pre- 
sented in this paper. However, experience with these transfer functions is still very 
limited, and care should be taken in extending their use to other particular applications. 

The results of the two examples studied show that the frequency response of the 
control system elements should be kept as high as possible to minimize the lag between 
the pilot's inputs and the system outputs. In the case of the lunar -landing simulator, it 
is concluded that the maneuvers performed with the simulator will involve dynamics 
which are possibly more difficult to control than the dynamics of the lunar module in the 
lunar environment. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 7, 1969. 
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APPENDIX A 


BRIDGE EQUATIONS 


The dynamic response of the hydraulic system was represented by the following 
equations: 


Qp - WpDp 


(Al) 


Qm _ w m^m 
Q L = ±k 0 a l |/p^| 

QrL = ±K o A Rv|/Pm| 


(A2) 

(A3) 

(A4) 


where Wp is determined by the dynamic equation 


T m _ K m ( w sync " w p) 


Thus 


and also 


Q c - Qp - Q m - Ql _ Qrl 


Q r = — p 


The relationships between hydraulic pressure and bridge velocity are then 




^m^m ~ Pm^m ^m^m n 


t l = Kt(^ m - i(, L ) 


R 


T L - B L 


(A5) 


(A6) 


(A7) 

(A8) 

(A9) 


Other nonlinear or dynamic factors included in the representation were (1) the time 
constant and limit displacement of the stroker and (2) selection of the gearing. 
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APPENDIX B 


CABLE MODES OF VIBRATION 

Sketch (b) shows the assumed representation for the cable, and the terms used in the 
derivation of the cable equations of motion: 



Sketch (b).- Cable and whiff letree. 


Assume that a solution exists in the form 


where 


y = h sin kat sin kz 


a 2 - _I_ 

w'/g 

T tension 


(Bl) 


w’ 


weight per foot (per meter) of cable 


A boundary equation, for the motion of the whiffletree, which must be satisfied is 


\dt y 


= -T 


z=V 


dz J z=i < 


+ ■ 


T y z=i' 


Manipulation of equations (Bl) and (B2) gives 

^_ mw k 2 a 2 + 3^t an k V = -Tk 


(B2) 


(B3) 
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APPENDIX B 


An iterative procedure was used to solve equation (B2) for k, where 

m w = 6.2 slugs (89 kg) 

w' = 0.8 lb/ft (1.2 kg/m) 

T = | m v = 8340 lb (3700 N) 

V = 200 feet (61 m) 

d = 6 feet (1.83 m) 

The value of k for the first mode of vibration was 0.015, and the frequency 
o> n = ka = 8.72 rad/sec. The Lagrange equation was then used to determine the coupling 
between the bridge motion, the pendulum motion, and the cable vibration. For the first 
mode of vibration, the equation for the motion of the cable and whiffletree is 

y = h sin 8.72t sin 0.015z 


and the velocity can be written 


y = h sin 0.01 5z 

the kinetic energy of the cable is, therefore, 


. , p200 , _ 

Kinetic energy = j j sin 0.015z + x + z 0 J‘ i d 2 


1.29h 2 + 3.28xh + 3440li + 2.48x 2 + 496x0 + 33 2OO0 2 


the kinetic energy of the whiffletree, located at z = 200, is 


Kinetic energy = ^ m w i[fi sin 0.015(200)J + x^ + 2000 > 


n 2 


= 0.06H 2 + 3. lx 2 + 12 4OOO0 2 + 0.86hx + 172h0 + 1240x0 


By applying the Lagrange equation, additional terms for the x and 0 equations, and 
the h equation were then determined 

|my + mg + 11.2^xg + (myZ + 1736^0 + 4.14h = K ^cp + K^< p (B4) 
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^rciyZ + 1736)x b + (myf2 + 314 4OO)0 + mygZ0 + 51 6h = 0 


(B5) 


2.70h + 4.14x + 5160 + 206h = 0 


(B6) 


where, according to Rayleigh’s principle, the coefficient of h in the h equation 
(eq. (B6)) has been set so that the uncoupled frequency of the first mode of vibration will 
be the same as that determined previously, 8.72 radians per second. 
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Figure 1.- Lunar-landing simulator. L-68-6959 

















deg 


h, ft 


*g, ft/sec 2 



- 2.3 


L / 

I I II I 


0.16 

0 


-i .04 

J- — o 


- 0.1 6 *— 

I — L_ I I I 


J-.04 



0 2 4 6 
Time, sec 


Figure 4.- Calculated cable divergence. 
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(a) Pilot-vehicle combination. 



(b) Pilot-vehicle-simulator combination with high gain simulator. 



(c) Pilot-vehicle-simulator combination with low gain simulator. 

Figure 6.- Calculated vehicle attitude response to a 225-foot (69 m) step translation command. 
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(a) Pilot-vehicle-simulator combination with low gain (4-second period) simulator. 
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(b) Pilot-vehicle-simulator combination with very low gain (8-second period) simulator. 

Figure 11.- Repeated calculation of vehicle attitude response to a 225-foot (69 m) step translation command. 
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(b) Pilot opinion rating as a function of w n f (from ref. 3, p. 57). 

Figure 12.- Effect on closed loop frequency and damping ratio of inserting first- and second-order servos in a single loop control task. 
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